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Two-particle jet-like azimuthal correlations have revealed intriguing modifications to the away- 
side of high p± trigger particles in relativistic heavy-ion collisions. Three-particle jet-like azimuthal 
correlation and three-particle azimuthal cumulant have been analyzed in experiments in attempt to 
distinguish conical emission of jet-correlated particles from other physics mechanisms. We investi- 
gate the difference between three-particle jet-like correlation and three-particle cumulant in azimuth. 
We show, under the circumstance where the away-side two-particle correlation is relatively flat in 
azimuth and similar in magnitude to the azimuthal average of the two-particle correlation signal, 
that the three-particle cumulant cannot distinguish conical emission from other physics mechanisms. 
The three-particle jet-like correlation, on the other hand, retains its discrimination power. 



PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 



Two-particle jet-like azimuthal correlations have re- 
vealed significant modification in central heavy-ion colli- 
sions at RHIC P, H i, i, i, i]. The away-side particles 
associated with and opposite to a high transverse mo- 
mentum (px) trigger particle are found to be broadly 
distributed about Acj) = tt in azimuth from the trig- 
ger particle, in contrast to observations from pp and 
d+Au collisions. The shape of the broad away-side dis- 
tribution varies with the associated particle p±. For 
1 < p± < 2 GeV/c, for instance, the away-side dis- 
tribution may even be double-humped with a dip at 
Acp = 7T @, S EL IE 13] • The away-side associated parti- 
cles are also found to be not much harder than the bulk 
medium particles [3, H, 0, [f|. The particles at Acj) — tt 
are found to be softer than those in the angular regions 
where the humps appear [3, El HI] , again in contrast to ob- 
servations in pp and d+Au collisions or jet fragmentation 
in vacuum. 

Several physical scenarios are possible to explain the 
observations. One is that jets may be deflected by radial 
transverse flow of the bulk medium or by the larger sur- 
vival probability of jet particles moving outwards than 
inwards due to energy loss Q. Such a scenario would 
have jet particles narrowly clustered in individual events 
but the cluster is randomly distributed about Acf> = tt 
over many events. The second is large angle gluon radi- 
ation Q. This scenario would have qualitatively similar 
structure as for deflected jets. The third is conical flow 
from sound shock-waves generated by large energy depo- 
sition of high momentum partons and its strong pressure 
disturbance in the medium [9(. Such shock- waves result 
in a distinctive Mach-cone type structure where particles 
are preferentially emitted at a Mach angle determined 
by the speed of sound in the medium, independent of the 
particle p± . If Mach-cone type conical flow is indeed re- 
sponsible for the observation, then the extraction of the 
speed of sound may be possible, thereby the equation 
of state of the created medium. The fourth is Cerenkov 



gluon radiation generated b y in teractions of fast moving 
particles with the medium [lCj. Such a scenario would 
have a similar structure as for the Mach-cone conical flow, 
but the Cerenkov angle will likely depend on the associ- 
ated particle p± [13] • 

Two-particle correlation cannot distinguish these sce- 
narios because they give qualitatively the same two- 
particle correlation. Such ambiguity is lifted in three- 
particle correlation and its p± dependence. If the broad 
two-particle correlation is due to deflected jets or large 
angle gluon radiation, the two associated particles will 
be narrowly clustered in angle but the cluster will swing 
over a wide range in azimuth on the away side. If the 
Mach-cone or Cerenkov radiation is responsible for the 
broad two-particle correlation, then the two associated 
particles will have equal probability being opposite away 
from Acj) = tt as being clustered together. A three- 
particle correlation signal with opposite azimuthal angles 
from A(j> = tt for the two associated particles is, there- 
fore, a distinctive signature of Mach-cone conical flow 
or Cerenkov radiation. The p± dependence of the cone 
angle may further discriminate between the two scenar- 
ios of Mach-cone conical flow and Cerenkov gluon radia- 
tion [13] . 

Three-particle jet-like azimuthal correlations have 
been studied in STAR 111] . The results show evidence 



of conical emission. The emission angle is independent 
of associated particle p± and is consistent with Mach- 
cone shock waves [3] , but not with simple Cerenkov gluon 
radiation [13]. The analysis followed the jet-correlation 
method commonly used in two-particle azimuthal corre- 
lation studies at RHIC [J [3, S El H, Q, but extended 
to three particles. Hereon we will refer to this method 
as the jet-correlation method. Recently another analy- 
sis method, the three-particle cumulant method with az- 
imuthal angle defined in the laboratory frame [13] , has 
been proposed. The method follows the mathematically 
well-defined cumulant concept. Hereon we will refer to 
this method as the lab-frame cumulant method. The two 
methods give different results and may confuse the gen- 
eral reader. In this paper we shall compare the two meth- 
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ods and discuss their differences in detail. We first give 
brief descriptions of the two methods. We then compare 
the two methods and discuss their differences, using a 
simple analytical model for jets. Finally we draw our 
conclusions. 



II. DESCRIPTIONS OF THE TWO ANALYSIS 
METHODS 

The objectives of the two analysis methods, the jet- 
correlation method and the cumulant method, are both 
to study jet structures. Due to the large particle multi- 
plicity in relativistic heavy-ion collisions, event-by-event 
reconstruction of jets is impossible; one often resorts to 
two- and three-particle azimuthal correlations of charged 
hadrons with high p± trigger particles that have a rel- 
atively large probability to originate from dijets. The 
obtained correlation functions, thus, yield information 
on dijets. In this section, we briefly describe the two 
analysis methods for three-particle azimuthal correlation 
studies. 



A. The jet-correlation method 

The three-particle jet-correlation method is described 
in [ll|, [ll|. The method is extended from the com- 
monly used two-particle jet-correlation method to three 
particles. Combinatorial backgrounds are obtained from 
event-mixing technique; they include background from 
three "random" particles as well as background from a 
correlated trigger-associated particle pair with a "ran- 
dom" associated particle. The "random" particles we 
refer to here (and hereafter without the quotes) may in- 
clude other correlations that are not related to the trigger 
particle, such as those due to anisotropic flow. 

The jet-like correlation method has the jet model 
in mind. The difficulty is that the underlying back- 
ground is unknown a priori. One has to make ad 
hoc working assumptions about the background level. 
The common assumptions made in data analysis are 
ZYA1 @, i, i, H, [II [3 and ZYAM The correla- 

tion measured at RHIC is the lowest around A(j> = ±1. 
The STAR experiment makes the assumption that the jet 
signal is zero within the fixed range of 0.8 < \A<fi\ < 1.2 
(ZYA1) 0, while the PHENIX experiment uses the so- 
called zero-yield-at-minimum (ZYAM) method in which 
the Acj) region where the signal minimum resides is de- 
termined by the data itself [4|]. 

The two-particle jet-like correlation is 

J 2 (A0) = J 2 (A0)-B 2 (A0), (I) 

where A(f> = <j) ~ <Pt is the azimuthal angle difference 
between associated and trigger particles. J 2 (A0) is the 
two-particle raw correlation function between the trig- 
ger and associated particles, and B 2 (A<fi) is the combi- 
natorial background normalized by the aforementioned 



normalization schemes. The three-particle jet-like corre- 
lation is 

J 3 (A0i,A</> 2 ) = J 3 (A0i,A^2)-B3(A^i,A</» 2 ) 

~J 2 (A^)B 2 (A^ 2 ) - J 2 (A0 2 )B 2 (A^), 

(2) 

or alternatively 

J 3 (A0 1; A0 2 ) = J 3 (A0 1 ,A^ 2 )- J B 3 (A^ 1 ,A0 2 ) 

-J 2 (A^ 1 )B 2 (A0 2 ) - J 2 (A<fe)B 2 (A0i) 
+2B 2 (Acj )1 )B 2 (A(j )2 ), 

where A<j)i = <f>i — <f>t(i — 1,2) are the azimuthal angles of 
the associated particles relative to that of the trigger par- 
ticle. J 3 (A^i, A<t> 2 ) is the three-particle raw correlation 
function. _B 3 (A0i , A(f> 2 ) is the combinatorial background 
of three-particle correlation between two random associ- 
ated particles with a random trigger particle. The third 
and fourth terms in the r.h.s. and Eqs. ([2]) and (J3]) are the 
other background, the combinatorial background of a cor- 
related trigger-associated pair with a random associated 
particle. It is given by the product of the two-particle 
jet-correlation signal J 2 with the underlying background 
B 2 . 

In data analysis 0, H, 0, H, @, EH , the dihadron back- 
ground B 2 (A(j)) is often obtained from the mixed-event 
technique, mixing a trigger particle from one event with 
an associated particle from another event, and then scal- 
ing the result by a normalization factor, a, using the 
ZYA1 or ZYAM scheme. The S 3 (A^i, A<fi 2 ) background 
is obtained from the mixed-event technique using three 
different events, with a proper normalization scale, a 2 b, 
where b quantifies the difference in the degrees of devia- 
tion from Poisson statistics in the background associated 
particle multiplicity and that from mixed-events (Tlj . In 
this paper, for simplicity we shall constrain ourselves to 
6=1. 

The advantage of the jet-correlation method is that, 
once the assumption about the background is made and 
the level of background is determined, the resultant 
three-particle jet-like correlation signal is easy to inter- 
pret and can be used to discriminate different physics 
scenarios. The disadvantage is of course the difficulty of 
the analysis. 

B. The lab- frame cumulant method 

The lab-frame three-particle cumulant method is de- 
scribed in detail in [l3[ . Given a high p± trigger particle 
to preferentially select a dijet, the two-particle cumulant 
is defined as 

P2(<k>4>) = P2 " Pl(<f>t)pi(<f>), (4) 

where <pt and <fi are azimuthal angles of the trigger and 
associated particles in the laboratory frame, respectively. 
The three-particle cumulant is defined as 

p3(<k><f>i,<h) = P3(<k><t>i,<h) - P2{<f>i,<h)pi(<k) 
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-p2(fa,fa)pi(fa) - p 2 (fa,fa)pl(fa) 

+2p 1 (fa) Pl (fa) Pl (fa), (5) 

where fa is the azimuthal angle of the zth associated 
particle (i — 1,2) in the laboratory frame. In Eq. |(SJ), 
Pi{4>) = dN/d(f> is the average single particle density, 
p 2 (fa,fa) = d 2 N / dfadfa is the average two-particle den- 
sity, and p 3 (fa,fa,fa) = d? N / dfad<j)\dfa is the average 
three-particle density. 

Cumulants are normally computed on a per event ba- 
sis. Since our goal here is to identify jet-correlation struc- 
ture, we normalize the cumulants by the number of trig- 
ger particles. The two normalizations differ simply by a 
constant factor which is the average number of trigger 
particles per event. 

It is worth to note here that pi, p2, and p3 are all 
event- wise average quantities. They depend on what 
event sample is used in the analysis. Consider for ex- 
ample two event samples within a given centrality bin: 
minimum bias event sample and only events containing a 
trigger particle (triggered events). The quantities related 
to a trigger particle, p 3 (fa, <j>i, fa), p 2 (fa,fa), p 2 (fa,fa), 
and p\(fa) are ah zero in the nontriggered events, and 
thus remain the same between the two event samples if 
normalized per trigger particle (or differ by a constant if 
normalized per event). The other quantities, p 2 (fa, fa), 
P i(fa), and pi (fa) are generally different between trig- 
gered events and nontriggered events due to trigger bias. 
Thus the cumulants analyzed using two different event 
samples differ, and this is beyond the simple normal- 
ization factor between per trigger and per event nor- 
malizations mentioned above. In other words, the cu- 
mulant analyzed using only the triggered events yields 
P3 (triggered events), the cumulant analyzed using only 
the nontriggered events would yield zero by definition, 
and the cumulant analyzed using both the triggered and 
nontriggered events together yields P3(all events), and 
P3 (triggered events) ^ p 3 (al\ events). Moreover, ps(all 
events) depends on the relative numbers of nontriggered 
events (for which p 3 = 0) and triggered events. In this pa- 
per, we shall constrain ourselves to using triggered events 
only. 

The per trigger normalized two-particle cumulant is 

p 2 (Afa = p 2 (Afa - Pl (Afa, (6) 

and the thrcc-particlc cumulant is 

p 3 (Afa,Afa) = p s (Afa,Afa)- p 2 (Afa,Afa) 

-p 2 {Afa) Pl (Afa) - p 2 {Afa) Pl {Afa) 
+2 Pl (Afa) Pl (Afa), (7) 

or alternatively 

p 3 (Afa,Afa) = p 3 (Afa,Afa)- p 2 (Afa,Afa) 

-p 2 (Afa) Pl (Afa) - p 2 {Afa) Pl (Afa). 

(8) 

In Eqs. ©, © and ©, p 2 (A0) = dN/dAcj) is the two- 
particle raw correlation function (or, equivalently, asso- 
ciated single-particle density per trigger particle in A<f> = 



(j) - fa), and p 3 (Afa,Afa) = (PN/dAfadAfa is the 
three-particle raw correlation function (or, equivalently, 
associated pair-density per trigger particle in A(f>i — 
fa - fa). In addition, p 2 (Afa,Afa) = p 2 (fa,fa) = 
d 2 N / dfadfa is the two-particle density but expressed in 
the azimuth differences of the particles from a random 
trigger particle. 

Because the azimuthal angles are defined in the labo- 
ratory frame, the single particle density p\(4>) is constant 
and we shall just use p\. 

pi(4>) = pi{Afa = pi = const. (9) 

The reason that pi(Afa) is also constant is because the 
product pi(fa)pi((j)) is taken after the pi's have been al- 
ready averaged over the event sample. 

Cumulants are well defined mathematically. The data 
analysis of cumulants is straightforward. One calculates 
the cumulants for each event and accumulates them over 
a event sample; the cumulants can be binned in fixed- 
multiplicity bins. The shortcoming, as we will see, is the 
difficulty in the interpretation of the obtained results: the 
cumulants are mathematically well defined and thus do 
not depend on the underlying physics of the events; their 
interpretations, therefore, have to depend on the model 
for the underlying physics. 

C. The reaction-plane-frame cumulant method 

Reaction plane direction in nucleus-nucleus collisions 
is random. The natural azimuthal angle to use is that 
defined relative to reaction plane, (j) — ip. The formulism 
in the previous section all applies with a simple change: 

<l)^fa-i) (io) 

and no change to A<p. However, P \ and p2 are not con- 
stant anymore, 

pi(4>) ^ const., 
pi(A0) ^ const., (11) 

but rather are functions of cf> — if) and A(f>, respectively. 

III. COMPARISONS BETWEEN THE TWO 
METHODS 

The lab-frame cumulant and the reaction-plane-framc 
cumulant are clearly different. This is because the cu- 
mulant is a mathematical construction, and an event 
sample expressed in lab-frame azimuthal angle and the 
same event sample expressed in azimuthal angle relative 
to reaction plane are different mathematical identities. 
However, the event sample contains obviously the same 
physics, whether expressed in the laboratory frame or 
in the reaction plane frame. In other words, the same 
physics yields two different cumulant results. Therefore, 
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extreme care should be taken in the physics interpreta- 
tions of cumulant results. 

The raw two- and three-particle correlation functions 
in the cumulant methods and in the jet-correlation 
method are identical, namely, 



p 2 (A0) = J 2 (A0), 



p 3 (A(f)i, A(f> 2 ) = J 3 (A0i, A&). 



(12) 



(13) 



Comparison of Eqs. ([3]) and ([7]) reveals that the reaction- 
plane-frame cumulant and the jet-correlation method are 
identical if the background B 2 (A(j)) and i? 3 (A0i, A(f> 2 ) 
are taken from mixed-events without normalization scal- 
ings (i.e., a = b = 1). This is because the same 
anisotropic flow is included in both methods. We shall 
therefore only focus on the comparison between the lab- 
frame cumulant method and the jet-correlation method. 

We shall first compare the two methods in a simple 
case, in which the background particle distribution is uni- 
form in azimuth (i.e., no anisotropic flow). Such a com- 
parison is enlightening as the difference between the two 
methods is straightforward and can be easily identified. 
We then include anisotropic flow in the comparison. We 
note that there could be intrinsic correlations (beyond 
those due to anisotropic flow) between the two associ- 
ated particles in the background. Those correlations are 
not important for our discussion here, so for simplicity 
we assume no other intrinsic correlations between the two 
associated particles except those from anisotropic flow. 



A. Jet-correlation with uniform background 

Since we consider no intrinsic correlation between the 
background particles and no anisotropic flow, we have 



p 2 (A(j>x, A0 2 ) = Pi{ 



i)=pl 



(14) 



Here we have assumed that the number of pairs equals to 
the square of the number of particles (Poisson statistics) . 
Equations ([7]) and ([5]) become 

/5 3 (A0i, A</) 2 ) = ps{A4>x, A^ 2 ) - pi{A<t>x)p\ 

-p 2 (Acf> 2 )p 1 +p 2 1 , (15) 

and 

/3 3 (A</>i, Acj> 2 ) = ps{A(j)x, A0 2 ) - p 2 {A4>x)px 

-p 2 {Afa)pi-p\. (16) 

Since no anisotropic flow is considered, the background 
distribution in jet-correlation method is uniform, 



B 2 (A<j ) ) = Bx, 



(17) 



where B\ is the average background single-particle den- 
sity. Since we consider no intrinsic correlation between 
the two background particles, we have 



Thus Eqs. © and © become 

J 3 (A^i,A0 2 ) = J 3 (A<t>x,A(f> 2 )-j 2 (A<Px)Bx 

-J 2 (A4> 2 )Bx-Bl (19) 

and 

J 3 (A0 1; A0 2 ) = J 3 (A<f>x,A(f> 2 )-J 2 (A<Px)Bx 

-J 2 {A4> 2 )Bx+Bl. (20) 

The single particle densities in the two methods are 
different: p\ is larger than B\, and 



Px-Bx = (J 2 ), 



(21) 



where ( j 2 ) is the average jet-correlated associated par- 
ticle multiplicity density in azimuth. As we shall see, 
this difference is the essential piece that makes the re- 
sults from the two methods differ. (Note that each event 
in our event sample contains a trigger particle. If events 
with no trigger particles are mixed into the event sample, 
then the difference px — Bx becomes rather arbitrary de- 
pending on the relative mixer and the specifics of trigger 
bias. See detailed discussion in the previous section.) 

Using Eqs. ©, dD, (HSJ), (HH), and (gT]), the difference 
between the three-particle cumulant, Eq. (I15[) . and the 
three-particle jet-correlation, Eq. (|2U)) . is given by 

A = p 3 (A<£i,A<fe)- J 3 (A</. 1 ,A0 2 ) 

= -(J 2 ) \j 2 {A<f>x) + J 2 (A0 2 ) - (J 2 )] . (22) 

To give further insights, we consider the simple case in 
which the three-particle jet-correlation function is simply 
the product of the two two-particle jet-correlation func- 
tions (i.e., the three-particle jet-correlation function is 
factorized) : 



J 3 (A^,A0 2 ) = J 2 {A4>x)J 2 (A<t> 2 ) 



(23) 



B 3 (A0i,A(fe) = fl?. 



(18) 



In other words, the two associated particles are not in- 
trinsically correlated, but are correlated due to their in- 
dividual correlations to the same trigger particle. This 
corresponds to the following extreme jet fragmentation 
scenario: the trigger particle direction is the jet axis, 
and fragmentations into individual hadrons are identical 
and independent of each other. Then, using Eq. (|2"2"1) . the 
three-particle cumulant is simply 

p 3 (A^,A0 2 ) = [j 2 (A0x) - (J 2 )] [J 2 (A0 2 ) - (J 2 ) 

(24) 

This should come as no surprise; the three-particle cu- 
mulant can also be factorized into two two-particle cu- 
mulants that are given by Eq. ©. The difference be- 
tween the two methods is in the background level. The 
jet-correlation method puts the background at ZYA1 or 
ZYAM, and the cumulant method, by definition, effec- 
tively takes the average of the raw correlation signal as 
background. In fact, if the jet-correlation method also 
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takes the average as the background, i.e., the two-particle 
jet-correlation signal is now J 2 (A0) — (J2) instead of 
J 2 (A</>), then it will yield the identical result as that from 
the cumulant method in Eq. (|24p . 

To give a visual comparison between the two methods, 



we use a specific example for the jet-correlation signal. 
We define a jet-like two-particle correlation with a near- 
side peak and a broad double-hump away-side distribu- 
tion: 



J 2 (A0) 



exp 



(A0) 2 
2a\ 



N 2 /2 



27TCT2 



exp 



(A(f> -7T + I 



2a 2 2 



\2l 



exp 



(A</> - 7T - (9) 5 

2^1 



(25) 



We suppose this jet-like two-particle correlation is 
present in every event (i.e., Mach-cone event), and it is 
atop a large uniform background, S true = 150/27T. We 
take the cone angle to be 9 = 1. We study two cases of 
jet signals: (A) the jets are narrow with o\ = 0.2 and 
(T 2 = 0.2; (B) we use a realistic jet signal as measured 
in experiment 0, H, [ll|, with o\ = 0.4 and er 2 = 0.7. 
For both cases, we take the numbers of jet-correlated 
particles to be TVi = 0.7 and N2 = 1.2, respectively, 
for near-side and away-side. These numbers are cho- 
sen so that the two-particle jet-correlation signal in the 
second case corresponds, qualitatively, to the measured 
one 0, [H, For both cases the minimum signal 

strength is at A<fi ss 1; for case (A) the normalized back- 
ground level is B± ss -Btrue (i-e., with a scaling factor of 
a = Bi/pi = 150/151.9), and for case (B) the normal- 
ized background is Si « -Btrue + 0.12 (i.e., with a scaling 
factor of a — B\j p\ — 150.12/151.9). For easy reference, 
we list the parameters below: 



(A): 



(B): 



Ni = 



N t = 

-Btrue 



0.7, N 2 = 1.2, (7i = 0.2, (j 2 = 0.2, = 1; 
= 150/27r,B 1 = B tIue , 



(26) 



0.7, N 2 = 1.2,(7! 

= 150/27T, Si = Strue 



0.4,(72 = 0.7,61 
0.12. 



The raw two-particle correlations are shown in Figs.QJa) 
and[2{a) for case (A) and (B), respectively, and will be 
described below in detail. 

We assume that the three-particle jet-like correlation 
can be factorized as the product of the two two-particle 
jet-like correlation functions given by Eq. (f2"5|) . The raw 
three-particle correlation is then given by 



r 2 

^truc 



J 3 (A^i,A^ 2 )= J 2 (A0 1 )J 2 (A0 2 ) 

+J 2 (A0i)S tme + J 2 (A<fe)S 

true 

J 2 (A^i) + S true ] [j 2 (A<fe) + S t 

(28) 

This raw three-particle correlation is needed as the start- 
ing point of any data analysis. One can extract the 
three-particle jet-correlation from the raw signal by mak- 
ing different assumptions of the underlying background 
level, which is the essential difference between the jet- 
correlation method and the cumulant method. 



The raw two-particle correlation for case (A) is shown 
in Fig. [DJa). The dotted line shows the real background 
level, Strue, (or the normalized background level Si, 
which equals to S tr ue)- The dash-dotted line shows the 
single particle density, p\ , used in the cumulant method. 
The three-particle correlation signal after subtraction of 
the real background (i.e., S trU e) is shown in Fig. QJb). 
This is the genuine three-particle correlation function 
that was initially put in, which shows the distinctive 
Mach-cone structure. The three-particle cumulant re- 
sult is shown in Fig. QJc) . This is equivalent to what the 
three-particle jet-correlation analysis would yield by sub- 
tracting an overestimated background level of Si = p\, 
given by Eq. flM} where ( J 2 ) = pi - B tluc . The Mach- 
cone structure is partially preserved in the three-particle 
cumulant result. The negative strips are produced by 
the over-subtraction of the average signal particle den- 
sity, which is larger than the true background level. 
More direct comparison between the two methods can 
be achieved by projections of the corresponding three- 
particle correlation signals. Figure Old) compares the 
projections along the off-diagonal axis on the away-side, 
A</>i ;2 > 1. The away-side on-diagonal projections are 
identical because the input to the simulation is symmet- 
ric. Both methods show a clear Mach-cone signal at the 
input cone angle (9 = 1). The structure from the cu- 
mulant result is more complex than that from the jet- 
correlation method. Our cumulant results are consistent 
with studies in [l3| where strong and narrow Mach-cone 
signal is also used. 

The raw two-particle correlation for case (B) is shown 
in Fig. [2fa). The dashed line shows the real back- 
ground level Btrue- The dotted line shows the nor- 
malized background level Si by ZYA1, such that the 
jet signal is zero at minimum. The dash-dotted line 
shows the single particle density, p\, used in the cu- 
mulant method. The three-particle jet-correlation sig- 
nal after subtraction of the real background level (i.e., 
Strue) is shown in Fig. [2jb). Again, this is the genuine 
three-particle jet-correlation function that was initially 
put in, J 3 (A0i,A(fe) = J 2 (A0i)J 2 (A0 2 ). Figure EJc) 
shows the three-particle cumulant result. The struc- 
ture of the cumulant result is quite complex and differ- 
ent from the three-particle jet-correlation result. This is 
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FIG. 1: (Color online) Uniform background case (A) with narrow jets and Mach-cones, Eq. (|26|l . (a) Two-particle correlation 
signal atop a uniform background. The dotted line is the estimated background level (Bi) to match the signal at the minimum 
(ZYA1), which is identical to the true background (Btrue)- The dash-dotted line is the average of the raw correlation signal 
(pi), (b) Three-particle jet-like correlation after subtracting the background (represented by the dotted horizontal line in the 
left panel), (c) Three-particle cumulant result, or three-particle jet-like result treating the average raw signal as background 
[represented by the dash-dotted line in panel (a)], (d) Comparison between the away-side off-diagonal projections of the three- 
particle jet-like correlation result in panel (b) and the cumulant result in panel (c). The on-diagonal projection is identical to 
the off-diagonal projection. The color bars on the right in panels (b) and (c) are three-particle correlation magnitudes. 
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FIG. 2: (Color online) Uniform background case (B) with realistic jets and Mach-cones, Eq. (|27[) . (a) Two-particle correlation 
signal atop a uniform background. The dotted line is the estimated background level (-Bi) to match the signal at the minimum 
(ZYA1), which is overestimated by SBi — 0.12 compared to the true background (-Btrue) plotted in the dashed line. The 
dash-dotted line is the average of the raw correlation signal (pi), (b) Three-particle jet-like correlation after subtracting the 
true background [represented by the dashed horizontal line in panel (a)], (c) Three-particle cumulant result, or three-particle 
jet-like result treating the average raw signal as background [represented by the dash-dotted line in panel (a)], (d) Comparison 
between the away-side off-diagonal projections of the three-particle jet-like correlation result in panel (b) and the cumulant 
result in panel (c). The on-diagonal projection is identical to the off-diagonal projection. The color bars on the right in panels 
(b) and (c) are three-particle correlation magnitudes. 



more clearly shown in the away-side off-diagonal projec- 
tions in Fig. [2^d) . The complex structure in cumulant is 
due to subtraction of the average single particle density, 
which is larger than the true background level. Unlike 
the simple case (A) where the jet peaks are narrow and 
well-confined, the broad jet peaks and over-subtraction 
of the background level in case (B) create the complex 
structure in the cumulant result. The distinctive Mach- 
cone structure on the away-side that is initially put in 
is hardly observable in the cumulant result. This is be- 
cause the away-side two-particle raw correlation as shown 
in Fig. HJa), qualitatively similar to that measured in 
real data [H, [f| [ll[ , is roughly flat and happens to have 
similar magnitude as the average single particle density. 
After subtraction of the single particle density, the away- 
side two-particle correlation strength is more or less sub- 
tracted away and can hardly show up in the final three- 
particle cumulant. 



The three-particle jet-correlation signal after subtrac- 
tion of the normalized background level (i.e., B\) is 
shown in the left panel of Fig. 02a). This would be 
the three-particle jet-correlation function from data anal- 
ysis using ZYA1 or ZYAM background normalization 
scheme. Since the background is overestimated, the ob- 
tained three-particle jet-correlation signal is lower than 
the true signal, but the structure of the genuine three- 
particle jet-correlation (i.e., the Mach-cone structure) is 
preserved. This is more clearly seen in the away-side off- 
diagonal projections shown in Fig.[3jb) where the dashed 
histogram is projected from the real three-particle corre- 
lation in Fig. [2tb) and the solid histogram is projected 
from the three-particle correlation result obtained with 
ZYA1 background normalization [i.e., from Fig. [3ja)] . 
In fact, the three-particle jet-correlation signal with the 
overestimated background, B\ > Struc, is given by 

J 3 (A0!,A0 2 ) = [j 2 (A0i) - 8Bi\ [i 2 (A^ 2 ) - 5B(h) 
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FIG. 3: (Color online) Uniform background case (B) with realistic jets and Mach-cones as in Fig. [5] (a) Three-particle jet-like 
correlation with background normalized (Bi = Btmc + 5B\, SB\ = 0.12) to the signal at A(f> = 1 [represented by the dotted 
horizontal line in Fig. HJa)]. Note the Mach-cone structure is still observable, (b) Away-side off-diagonal projections of the 
three-particle jet-like correlation result subtracting the true (input) background [shown in Fig.^b)] and that subtracting ZYA1 
background [shown in panel (a)]. The dashed curve is as same as that in Fig. [2jd). (c) Zoomed-in version of Fig. [2jc). The 
color bars on the right in panels (a) and (c) are three-particle correlation magnitudes. 



(where 5B\ = B%- B truo ). 



Note the similarity between Eqs. (j2"5|) and (|2"4"j) : if the uni- 
form background levels are chosen to be equal, then the 
resultant three-particle jet-correlation and three-particle 
cumulant are identical. 

On the other hand, if one zooms into the cumulant 
result in Fig.[^c), now shown in Fig. [3Jc), one may also 
see the Mach-cone structure due to the double-hump in 
the away-side two-particle correlation. However, it would 
be extremely difficult to distinguish it from the many 
other and much larger peaks in the cumulant result. 



B. Anisotropic flow without jet-correlation 



where and v 2 are the elliptic flow parameters of the 
trigger and associated particles, respectively. Likewise 
i?4 and «4 are the respective fourth harmonic coeffi- 
cients. The two-particle jet-correlation signal of Eq. (Q]) 
is then, 



J 2 (A0)= J 2 (A0) 
-B 1 



2v { 2 ] v 2 cos(2A<?(>) + 2v[ t) v 4 cos(4A0) 

(31) ' 



In heavy-ion collisions, particles are correlated to re- 
action plane due to the hydrodynamic type of collective 
flow of the bulk medium and the anisotropic overlap re- 
gion between the colliding nuclei. The trigger particle 
emission is also correlated to reaction plane due to, not 
so much of hydrodynamic flow, but the path-length de- 
pendent energy loss of high p± particles in the medium 
that is initially anisotropic. This reaction plane correla- 
tion, expressed in harmonics up to the fourth order, is 
given by 



B 2 {Acj)) = B x 1 + 2v [ 2 ] v 2 cos(2A0) + 2v 



V4 cos(4A</>) 
(30) 



There are two combinatorial backgrounds to three- 
particle jet-correlation. One is that of a correlated 
trigger-associated pair combined with a random back- 
ground particle. This background can be obtained 
by the product of the two-particle jet-correlation sig- 
nal with the underlying background particle, namely 
J 2 (A</>i)B 2 (A0 2 ) + J 2 (A^ 2 )S 2 (A0i). The other combi- 
natorial background is due to the trigger particle com- 
bined with two random particles from the underlying 
background. This background term normalized per trig- 
ger particle, considering only the anisotropic flow corre- 
lation, is given by 12] 



B a (A&,Afo) = (Bj 



( 1 +2v ( 2 t) v ( 2 1) cos(2A^i) + 2v { 2 t) v { 2 ) cos(2A0 2 ) + 2v i 2 1) v { 2 ) cos2(A0i - A0 2 ) \ 

+ 2v [ l ) v [ l ) cos4(A0! - A0 2 ) 



-2v 
-2v. 



V +2v 



cos(4A0! 



2vfvf ] cos(4A0 2 ) 



cos2(A0i 
' cos 2(A(/)i 



2A0a) 
A0 2 ) 
I 



2v 



(t)„,(2)„,(l) 



cos2(2A0i - Acf> 2 



(32) 



where (Bf) is the per trigger associated pair density which is not necessarily equal to (Si) 2 = B\ (i.e., non- 
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FIG. 4: (Color online) Lab-frame three-particle cumulant from pure anisotropic flow correlation, Eq. (|36|l , (a) Non-Poisson 
statistics induced cumulant with typical (Si) — (Bi) 2 = O.l(Bi). (b) Irreducible flow correlation terms (where (B 2 ) — (Bi) 2 is 
taken), (c) Total cumulant given by Eq.|35][sum of panel (a) and (b)]. (d) Away-side on-diagonal and off-diagonal projections of 
the total cumulant in panel (c). The color bars on the right in panels (a), (b), and (c) are three-particle correlation magnitudes. 



Poisson statistics). In Eq. (132|) and hereafter, and 

(2) 

v 2 are elliptic flow parameters of the two associated 

particles, and and ti^ are the respective fourth har- 
monic coefficients. The final three-particle jet-correlation 
function is given by Eq. ([2]) or Eq. (J3]) . 

It is interesting to consider only anisotropic flow cor- 
relation, no jet-like correlation. Obviously, if no jet-like 
correlation is present, the three-particle jet-correlation 
will give zero signal as there will be no two-particle jet- 
correlation signal to start with, J 2 (A</>) = 0. However, 
the three-particle cumulant will still yield nonzero result, 
as we demonstrate below. 

With only anisotropic flow correlation present, the raw 



one-, two-, and three-particle lab- frame cumulants are 
simply given by the background, p\ = B±, p 2 {A<p) = 
B 2 (A(j)), and p 3 (A0i, A<j) 2 ) = B 3 (A0i, A<fe). Thus the 
three-particle cumulant from Eq. ([7]) is 



p 3 (A0!, A0 2 ) = B 3 (A0i, A0 2 ) - B 2 (Acj) 1 - A<fe) 

-B 2 (A ( j )1 )B 1 - B 2 (Acj) 2 )B 1 + 2B 2 . 

(33) 

Here B 2 {A4>) is given by Eq. ((301), B 3 (Acp 1 , A(j) 2 ) by 
Eq. (f3"2")> . and B 2 (Acf>i — Acf> 2 ) is the two-particle flow 
correlation normalized per trigger particle, and is given 
by 



B 2 (Afa - A0 2 ) = (B\ 
With simple algebra, we obtain 

I O/Dl 



- 2v ( 2 ] v ( 2 ] cos2(A^i - A4> 2 ) + 2v { l ] v [ p cos4(A0! - A0 2 )j . (34) 



+2{Bf) [4 t3 4 13 4 2) 



4 t} ^ 1} cos(2A0!) + u^u^ cos(4A0i) + v { 2 ] v { 2 ] cos(2A0 ; 
! -n2(A<£i - 2A</> 2 ) + i^fM* 3 cos2(2A0i - Ac/> 2 ) + ? 



vf>vf> cos(4A</> 2 



'4 cos 



> 4 1} cos2(2A0i - Afc) + v^v^vf cos2(A0i + A<j) 2 ) 
(351 



With no jet correlation, only anisotropic flow correla- 
tion present, the three-particle cumulant gives non-zero 
correlation result. The three-particle lab- frame cumulant 
measures something different from simple jet-correlation; 
it measures three-particle correlation regardless of the 
nature of the underlying physics. With Poisson statis- 
tics, (B 2 ) = (Bi) 2 , the first square-bracket term in the 
r.h.s. of Eq. (|35() drops and what remains are the ir- 
reducible flow correlation terms (in the second square 
bracket). Although these terms are on the order of v 2 , 
their magnitudes can be still sizable due to the large 
background level of B\ in relativistic heavy-ion collisions. 
With non-Poisson statistics, the cumulant is even more 
complicated. To give a visual impression of the cumulant 



from pure anisotropic flow, we plot Eq. (|35[) with typical 
flow magnitudes: 

vf = 7.5%, v { 2 1] = 4 2) = 5%,u 4 = vl (36) 

The result is shown in Fig. 2] Panel (a) shows the non- 
Poisson statistics induced cumulant with typical 

(B?) - (Bi) 2 = 0.1(5!). (37) 

Panel (b) shows the irreducible flow correlation terms 
where we have taken (B 2 ) = (Bi) 2 . Panel (c) shows 
the total cumulant. Panel (d) shows the away-side on- 
diagonal and off-diagonal projections of the total cumu- 
lant. As seen, pure flow gives non-zero cumulant with 
complex structures. 
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C. Jet-correlation with anisotropic flow 
background 

We now turn our attention to the situation where both 
jet-correlation and anisotropic flow are present. For the 
sake of simplicity we shall assume Poisson statistics 

(Bl) = (B x ) 2 = Bl (38) 

If the true background level and the anisotropic flow 
magnitude are precisely known, then all the background 
terms are determined, and the final jet-like three-particle 
correlation signal from the jet-correlation method recov- 
ers the true input signal of Eq. ([25| . the same as that 
shown in Fig. [2(b). However, experimentally the true 



background level is not known a priori, and thus back- 
ground normalization schemes, such as ZYAI or ZYAM, 
have to be employed. We shall examine the effect of 
an imperfect background estimate in the jet-correlation 
method. We will not discuss the effect of uncertainties 
in the anisotropic flow magnitudes as they have been ex- 
tensively discussed in Ref. 12]. Moreover, systematic 
uncertainties on the flow measurement can be reliably 
accessed in experiment. 

From Eqs. ©, (O, (EH), and (02), an uncertainty 5B X 
in the background level B\ introduces a change in the 
three-particle correlation signal from the jet-correlation 
method: 



SJ 3 



[2B 1 8B 1 + (SB!) 2 } 



( 1 +2v ( 2 ] v ( 2 ] cos(2A0i) + 2v { 2 ] v ( 2 2) cos(2A<?(> 2 ) + 2v ( 2 1) v ( 2 2) cos2(A0i - A</> 2 ) \ 
+2vf ) v ( 4 1) cos(4A0 1 ) + 2vf ) vf ) cos(4A0 2 ) + 2v ( 4 1) v (2) cos4(A0! - A</> 2 ) 



2A</> 2 ) + 2v K 2 'v. 



+2v^ , v ( 2 L, vX > cos2(A0i 
V +2v { 2 1] v {2) v [ 4 t] cos 2(A0j + A0 2 ) 

-(SBi J 2 (A0i) [l + 2v$ ) vf ) cos(2A<£ 2 ) + 2v<£ ) vf ) cos(4A</> 2 ) 



(t)„,(2) fl) 



cos2(2A(/>i - A(f) 2 ) 



5BiJ 2 (A0 2 ) 1 + 2v [ 2 ) v ( 2 1) cos(2A(/. 1 ) + 2v%>v\ L> cos(4A^ 1 



2[B 1 SB 1 + [5Bi) z ] 1 + 2v { 2 ] v {2) cos(2A^ 2 ) + 2v% J v£ J cos(4A0 2 ) 1 + 2v) 2 'v y 2 > cos(2A0i) + 2v^vX' cos(4A^) 



.,(*)., ( 2 ) 



/ 



-5B 1 J 2 (A(j) 1 ) 1 + 2v { 2 ] v {2) cos(2A0 2 ) + 2v% , vf cos(4A<£ 2 ) 
-5B 1 J 2 (A<h) \l + 2v i 2 ) v { 2 1) cos(2A0i) + 2v i 4 t) v i 4 1) cos(4A0i) 



(t) (2) 

\ / „, \ I 



-(SB,) 2 



-\2Bx5Bx 



1 + 2v [ 2 ) v < 2 2> cos(2A0 2 ) + 2v 4 z, v 4 z > cos(4A0 2 



2v ( 2 t) v { 2 1) cos(2A<?!)i) + 2v£>v1 ) cos(4A0i) 



/ 2v. 



(!)„,(2) 



cos2(A0i - A</> 2 ) + 2v[ 1) v {2) cos4(A0i - A^ 



(^i) 2 ] 



\ 



2v£>v { 2 1) vf 5 cos2(A</>i - 2A</> 2 ) + 2u 2 t) u 2 2) 4 1) cos2(2A0i - A<£ 2 )+ 
2v [ 2 1) vf ) v { * ) cos2(A<fc + Acf> 2 )- 

4(u 2 t) ) 2 4 1) w 2 2) cos(2A0i) cos(2A0 2 ) - Av { 2 t] v ( 4 ] v ( 4 1] v ( 2 2) cos(4A</>i) cos(2A</> 2 )~ 
V 4i4 t M t) 4 1 M 2) cos(2A&) cos(4A0 2 ) - 4(4 t) ) 2 ^ 1) vf ) cos(4A<£i) cos(4A0 2 ) / 



(39) 



If the true three-particle jet-correlation signal is given by Eq. 
ground normalization scheme) becomes 



, then that with the scaled background (by back- 



J 3 (A</>i,A0 2 ) 



J 2 (A0i) 
J 2 (A0 2 ) 



5Bi 
8Bi 



1 + 2v { 2 ] v { 2 ] cos(2A0i) + 2v ( l ] v { l ] cos(4A0i) 
1 + 2v i 2 ) v {2} cos(2A0 2 ) + 2v ( A t] v {2) cos(4A<?(> 2 ; 



-\2Bx5Bx + (SBxf 



( 2v [ 2 ) v { 2 ) cos2(A</>i - A0 2 ) + 2v ( i 1) v (2) cos4(A<?(>i - A</> 2 )+ \ 
2?4 t) « 2 1) «4 2) cos2(A0i - 2A0 2 ) + 2v ( 2 t) vf 2 2) v[ 1) cos2(2A(^ - A0 2 )+ 
2^ 1) 4 2) «i t) cos2(A0x + A0 2 )- (4 3) 

4(^ t) ) 2 ^ 1) ^ 2) cos(2A<?(>i) cos(2A0 2 ) - 44* ) ^ t) 4 1) ^ 2) cos(4A0 x ) cos(2A0 2 )- 
V 44 t) 4 t) 4 1) ^4 2) cos(2A0!) cos(4A0 2 ) - 4(^ t) ) 2 4 1) vf ) cos(4A<?!> 1 ) cos(4A0 2 ) / 



We postpone the discussion of the above result to later. We now proceed to the calculation of the lab-frame 
three-particle cumulant. From Eqs. 0, (TT2")) . (TIB"]) , (j2Tj) , and using 

p 2 (A0i, A0 2 ) = p\ [l + 2w 2 cos2(A& - A0 2 ) + 2w 2 cos4(A0! - A0 2 )] = ( / o 1 /B 1 ) 2 B 2 (A< ? !) 1 - A<fe), (41) 
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we obtain 

p 3 (A0i, A0 2 ) = J 3 (A0i, A02)-[J 2 (A0i) + J 2 (A0 2 )] (s x + (J 2 ))-(l + <J 2 )/B 1 ) 2 S 2 (A0 1 -A0 2 )+2 (s x + <J 2 )) 
Taking the difference between Eqs. and ©, using Eqs. ([521 and (p?4"|) , we have 



(42) 



A = ^(A^A^)- J 3 (A0i,A<fe) 
= J 2 (A^i)(£i 



(A^ 2 ) ( 



Bi 



1 + 24 t) «2 2) cos(2A0 2 ) + 2vfvf ) cos(4A</> 2 ) - B 1 - (J 2 ) 
1 + 2«|* ) 4 1) cos(2A^i) + 24 t) «i 1) cos(4A0i)l - Si - (J 2 ) 



1 +2< ; 4 cos(2A0i) + 2v\ t} v 2 > cos(2A0 2 ) + 2u 2 j 4 cos2(A0i - A0 2 ) 
+2v { i ) v [ l ) cos(4A0i) + 2^ t} 4 2) cos(4A0 2 ) + 2v { pv ( p cos4(A0i - A0 2 ) 

+2v < £ ) v < 2 V) vf ) cos2(A< ? !) 1 - 2A0 2 ) + 2v { £ ) vf ) v ( £ ) cos2(2A</>i - A0 2 ) + 2v 2 1) v < £ ) vf ) cos2(A0i + A</> 2 ) 



- (B'f + 2B X {J 2 ) + {J 2 ) 2 ) 1 + 2v < £ ) v < £ ) cos2(A<^ - A<£ 2 ) + 2v^>v^' cos4(A0! - A0 2 ) 



(!)„,( 2 ). 



+2 (^ J B 1 2 + 2 J B 1 (J 2 ) + (J 2 ) 2 

-2B 2 [l + 2v [ 2 t) v ( 2 1) cos(2A0i) + 24%^ cos(4A0i) 1 + 2v ( 2 t) vf > cos(2A</> 2 ) + 2v ( f ) vf > cos(4A0 2 ) 
= -<J 2 ) (j 2 (A0 1 ) + J 2 (A0 2 )- (J 2 )) 

(i 2 (A0x) - (J 2 )) [2«^4 2) cos(2A0 2 ) + 2v { i ) vf ) cos(4A<fe) 
+Bi (j 2 (A0 2 ) - (J 2 >) [2^ t) 4 1) cos(2A0 1 ) + 2v[ t) v { 4 1) cos(4A0 1 ) 
-<J 2 ) (2Bi + <J 2 )) [24 1} 4 2) cos2(A0! - A<fe) + 2v[ 1) v {2) cos4(A0! - A<fe) 

+B\ hv^v^vf) cos2(A0i - 2A</> 2 ) + 2v ( 2 ] v {2) cos2(2A0i - A0 2 ) + 2v { 2 1) v { 2 2) v[ t] cos2(A0i + A0 2 )j . (43) 
Assuming factorization of the three-particle correlation signal in Eq. (|25p , the three-particle cumulant is therefore 



p 3 (A0i,A0 2 )= J 2 (A0 1 )-(J 2 ) J 2 (A0 2 )-(J 2 ) 



+B 1 t Ja(A^i) - (J 2 ) 
+Bi (J 2 (A0 2 ) - (J 2 ) 



2t4 t) *4 2) cos(2A</> 2 ) + 2^?;^ cos(4A<?(> 2 ) 
2v ( 2 ) v ( 2 1) cos(2A0i) + 24 t) 4 1) cos(4A<^) 



(t) a ,(2) 



-(J 2 ) f 2fli + (J 2 )J cos2(A0! - A0 2 ) + 2< J ^ J cos4(A0! - A0 2 

+B\ hv^ v { 2 } v {2) cos2(A0i - 2A0 2 ) + 2v { 2 } v {2) v[ 1} cos2(2A0i - A</> 2 ) + 2v { 2 ) v (2) v[ t] cos2(A0i + A0 2 ) 

(Ml 



(!)„,( 2 ) 



The hrst line in the r.h.s. of Eq. (|44|) is the three-particle 
cumulant result for the case of no anisotropic flow cor- 
relation, and the fifth line is the three-particle cumu- 
lant due to irreducible anisotropic flow correlation [as in 
Eq. 1(55])] . However, there are additional terms in the rest 
of Eq. (|44]) r.h.s., which arise from the coupling between 
jet-correlation and anisotropic flow. 



We are now ready to compare the jet-correlation 
method and the lab-frame cumulant method for the real- 
istic case of jet-correlation together with anisotropic flow 
background. We again use the jet-model of case (B) in 
Eq. (|27[) . but with additional anisotropic flow of magni- 



tudes given by Eq. 



Namely, 



Ni = 0.7, N 2 = 1.2, o-i = 0.4, cr 2 = 0.7, ( 
B x = Struc = 150/27T, SBx = 0.12; 
v { 2 t] = 7.5%, ^ = v (2) = 5%, Vi = vl 



(45) 



We first examine the effect of an imperfect background 
estimation (by ZYAl or ZYAM normalization) on the fi- 
nal three-particle jet-correlation signal given by Eq. (|40|) . 
Figure EJa) shows in the solid curve the raw two-particle 
correlation, J 2 (A0) given by Eqs. (|3"Tj) and (|2"5")) with pa- 
rameters in Eq. (|4"5"|) , and in the dashed curve the flow 
modulated background of Eq. (f30|) . The dotted curve 
shows the scaled background by ZYAl. Figure Ob) 
shows the final three-particle jet-like correlation obtained 
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FIG. 5: (Color online) Jet-correlation model with realistic jets and Mach-cones together with anisotropic flow background, 
Eq. (|45(l . (a) Two-particle correlation signal, J2(A0) by Eq. (I31|l (solid curve) atop a flow modulated background, £>2(A<^) 
by Eq. (|30[) (dashed curve). The dotted curve is the (over)estimated background (by 8B1 = 0.12) to match the signal at the 
minimum (ZYA1). The dash-dotted line is the average of the raw correlation signal, pi. (b) The three-particle jet-correlation 
after subtraction of the ZYAl-normalized background [represented by the dotted curve in panel (a)], (c) The lab-frame three- 
particle cumulant. (d) Away-side on- and off-diagonal projections of the jet-like result in panel (b) in thin curves and cumulant 
result in panel (c) in thick curves. The color bars on the right in panels (b) and (c) are three-particle correlation magnitudes. 
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FIG. 6: (Color online) Jet-correlation method details using jet-correlation model with realistic jets and Mach-cones together 
with anisotropic flow background, Eq. (I45|l . (a) The error introduced by an inaccurate background level, the second term 
in r.h.s. of Eq. (|40[) . (b) Illustration of what the constructed three-particle correlation would be if the over-subtracted jet- 
correlated particles have the same anisotropic flow as the medium background particles so that the inaccurate background 
would not introduce structures shown in the left panel. The three-particle jet-correlation would be similar to the input in 
shape with a reduced magnitude, and would be similar to that constructed in the uniform background example shown in 
Fig. [2jb). (c),(d) Away-side off-diagonal (c) and on-diagonal (d) projections of the input Mach-cone signal [i.e., as same as 
that in Fig. [2Jb)] in dash-dotted curves, of the structure of the flow effect in panel (a) in dotted curves introduced by an 
over-subtracted background level, of the three-particle signal in panel (b) in dashed curves that would be constructed if no 
structure was introduced by the inaccurate background level, and in solid curves of the obtained final three-particle signal using 
ZYA1 background subtraction in Fig. [SJb). The solid curves in (c) and (d) are as same as the thin dashed and solid curve in 
Fig.EJd), respectively. The color bars on the right in panels (a) and (b) are three-particle correlation magnitudes. 



by Eq. after subtraction of the scaled backgrounds. 
The result is significantly distorted from the input cor- 
relation [as shown in Fig.[5fb)], however, the Mach-cone 
structure seems still evident. The distortion is due to 
flow effect caused by subtraction of an incorrect back- 
ground level by ZYAl. Figure [5jd) shows the away-side 
on- and off-diagonal projections in thin curves. The on- 
and off-diagonal projections are not the same any more; 
the symmetry in the input Mach-cone signal is lost. The 
flow effect due to over-subtraction of the background level 
is different between on- and off-diagonal projection; the 
shape of the on-diagonal projection is not affected but 
that of the off-diagonal project is (also see below). The 
net effect is that the peaks in the off-diagonal projection 
are pulled toward the flow peaks at ±7r/4. The effect of 
flow due to incorrect subtraction of the background level 
will be further discussed below. 



Figure [HJc) shows the lab-frame three-particle cumu- 
lant by Eq. (|44]) . The cumulant result is very differ- 
ent from the three-particle jet-correlation result shown 
in Fig. [2Kb) or Fig. \5\b). The structure of the cumu- 
lant is complex. Figure E{d) also shows the away-side 
on-diagonal and off-diagonal projections of the cumulant 
result in thick curves. The on- and off-diagonal projec- 
tions are out of phase; the on-diagonal projection is even 
all negative. The input Mach-cone signal, which consists 
of two off-center peaks in both on- and off-diagonal di- 
rections, is not observable. 

The distortion of the constructed jet-like three-particle 
correlation is due to the overestimated background level 
by 8B1/B1 ~ 0.5%. The jet-correlation signal from 
Eq. (|25|) does not contain anisotropic flow. However, 
part of the signal is now included into the background 
which is all taken to be flow-modulated. This introduces 
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flow structures in the over-subtracted part of the three- 
particle background which should be really part of the 
signal that is uniform. The reverse of the over-subtracted 
background is given by the last line of r.h.s. of Eq. (|4"0")) ; 
we show this in Fig.JSJa). 

If the number of jet-correlated particles also varies with 
reaction plane in the same way as the medium particles, 
then the constructed three-particle jet-correlation would 
be similar to the input correlation except an approxi- 
mately constant reduction in the signal magnitude. In 
other words, the final signal would be as same as given 
by the first line of r.h.s. of Eq. (|40|) . This is shown in 
Fig. Hfb). The signal shape is not significantly altered 
from the true signal in Fig. [2fb). The argument can be 
turned around: if the medium particle background is uni- 
form relative to reaction plane, then the over-subtracted 
background due to ZYA1 normalization will not intro- 
duce structures and the final three-particle correlation 
would be similar to the true one except a constant re- 
duction in the signal magnitude [Fig. Ela)]. The results 
in Figs.[3fa) and[6^b) are indeed very similar. 

Figures [D(c) and (d) show the effect of flow due to 
subtraction of an incorrect background level in off- and 
on-diagonal projection, respectively. The dash-dotted 
curve shows the input Mach-cone signal [projection of 
Fig. 1211b)]; the dashed curve shows the projection of 
the jet-like three-particle correlation by ZYA1 back- 
ground subtraction but excluding flow effect [projection 
of Fig. [U[b)] ; the dotted curve shows the projection of the 
flow effect in Fig^a); and the thick solid curve shows 
the projection of the final result [projection of Fig. [5jb)] , 
which is the sum of the dashed and dotted curves. The 
ZYA1 background subtraction causes an reduction in 
the three-particle correlation magnitude. The flow effect 
changes the magnitude of the on-diagonal projection, but 
does not affect its shape. However, the flow mismatch be- 
tween jet-correlated particles (no flow in our particular 
model example) and medium background particles (finite 
flow in our model example) affects off-diagonal projection 
in both its magnitude and shape. The flow effect pulls 
the off-diagonal peak positions from the input ±1 toward 
the flow peaks at ±rr/4. 

In reality, the situation may be somewhere in-between: 
the jet-particles do have some variation relative to reac- 
tion plane, and the variation may not be as strong as 
that of the medium particles. If this is true, then the 
constructed three-particle jet-correlation signal should 
be somewhere in-between those in Figs. EJb) and[6jb). 
We try to illustrate this by subtracting half of the jet- 
correlation pedestal as background (i.e., setting SBi = 
0.06), effectively modulating the jet-correlation pedestal 
with half of the medium background flow. The result is 
shown in Fig.JTJa). The Mach-cone signal is less distorted 
than in Fig.^b), and the cone angle is closer to the in- 
put 8=1. However, if jet-particles are more strongly 
modulated than background particles (due, for instance, 
to the recently observed ridge particles then the 

over-subtracted flow modulation by ZYAl normalization 



would be too weak. The resulting jet-correlation shape 
may be mimicked in our jet-correlation model by sub- 
tracting a lower level of background (i.e., 5B\ < 0). 
We illustrate this in Fig. [TJb) using 5B 1 = -0.06. As 
expected, the off-diagonal peaks are pushed away from 
the input angles ±1 toward ±7r/2 (instead of ±7r/4 for 
5Bi > 0), while the on-diagonal peak positions remain 
intact. On the other hand, it is also possible that the 
jet-particle modulation with respect to reaction plane is 
opposite to that of background particles, i.e., jet-particles 
have a negative elliptic flow parameter. This can come 
about at relatively low p± due to the path-length effect 
of jet- quenching; more low p± particles are generated in 
the direction perpendicular to reaction plane where more 
high p± particles are quenched. Such a scenario would 
make the flow effect due to ZYAl background normaliza- 
tion more significant than shown in Fig. [5](b) , resulting 
in off-diagonal peaks closer to the flow peaks at ±7r/4. 

If jet-correlation signal does not vary with reaction 
plane angle (as in our jet-correlation model), then the 
evolution from Fig. [TJb) (5B X = -0.06) to Fig. Hfb) 
(SBx = 0) to Fig. (Ha) (SB-l = 0.06) to Fig. EJb) 
(6Bi = 0.12) gives a good indication of the background 
normalization effect. The magnitude drops because of 
the larger background subtraction, and the shape evolves 
because of the flow effect due to subtraction of incor- 
rect background. This evolution is better recapitulated 
in Figs. [71c) and (d) which display, respectively, the off- 
and on-diagonal projections of the jet-like three-particle 
correlation signals obtained with different levels of back- 
ground subtraction. The increasing flow effect with in- 
creasing SB\ pulls the off-diagonal peaks from the input 
cone angle ±1 toward the flow peaks at ±tt/4, while an 
under-subtraction of background (SBi < 0) pushes the 
off-diagonal peaks away towards ±7r/2. The on-diagonal 
projection, on the other hand, only changes its magni- 
tude but not its shape because the flow effect is unifrom 
along the on-diagonal direction as indicated in Figs. UK a) 
and (d). 

Experimentally, the extracted conical emission angle 
is found not to be strongly affected by uncertainties in 
the background normalization level This suggests 

that the jet-correlation signal has variations with reaction 
plane similar to the medium background. 

The complex structure of the lab-frame three-particle 
cumulant shown in Fig.[5Jc) is due to the irreducible flow 
correlation terms and the jet-flow cross-terms. The irre- 
ducible flow correlation terms, Eq. (|35[) . is already shown 
in Fig.HJb). The contribution of the cross-terms between 
jet-correlation and anisotropic flow is shown in Fig.[5Ja). 
The three-particle cumulant excluding the jet-flow cross- 
terms is shown in Fig. [U(b) . Likewise, the three-particle 
cumulant excluding the irreducible flow correlation terms 
is shown in Fig. |5fc). The three-particle cumulant ex- 
cluding both the irreducible flow terms and the jet-flow 
cross-terms is show in Fig. [Hd). The corresponding 
on- and off-diagonal projections are shown below each 
two-dimensional plot in panels (e)-(h), respectively. Of 
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FIG. 7: (Color online) Jet-correlation method background normalization effect using jet-correlation model with realistic jets 
and Mach-cones together with anisotropic flow background, Eq. (|45[) . (a) The three-particle jet-like correlation signal after 
over-subtraction of a background half way to ZYA1 (i.e., by SBi = 0.06). (b) The three-particle jet-like correlation signal after 
under-subtraction of background by SBi = —0.06). (c),(d) Away-side off-diagonal (c) and on-diagonal (d) projections of the 
jet-like correlation signals using different background normalizations by SBi. The thick dashed curves (SBi — 0) are the input 
signal [and the same as the dashed curve in Fig. Eld)], and the thick solid curves (SBi — 0.12) are the final signal by ZYA1 
normalization [and the same as the respective solid curves in Fig. [6jc,d)] . The color bars on the right in panels (a) and (b) are 
three-particle correlation magnitudes. 
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FIG. 8: (Color online) Cumulant method details using jet-correlation model with realistic jets and Mach-cones together with 
anisotropic flow background, Eq. (|45[l . (a) The jet-flow cross-terms in lab-frame three-particle cumulant [i.e., sum of lines 2, 3, 
and 4 of the r.h.s. of Eq. (|44|l ]. (b) The lab-frame three-particle cumulant excluding the jet-flow cross terms [i.e., sum of lines 1 
and 5 of the r.h.s. of Eq. (|44|l]. (c) The three-particle cumulant excluding the irreducible pure flow correlation terms [i.e., sum 
of lines 1, 2, 3, and 4 of the r.h.s. of Eq. (|44[) ] . (d) The lab-frame three-particle cumulant excluding both the jet-flow cross-terms 
and the irreducible flow terms [i.e., line 1 of the r.h.s. of Eq. (144 [) ]. This is similar to that constructed in the uniform background 
example shown in Fig. [5Jc). (e)-(h) Away-side on- and off-diagonal projections of panels (a)-(d), respectively. The projections 
are identical in (h). The color bars on the right in panels (a), (b), (c) and (d) are three-particle correlation magnitudes. 



course, what is shown in Figs. E^d) and (h) is now in 
the same spirit of the jet-correlation method, subtract- 
ing specific background terms. Even with that, the input 
Mach-cone signal is not observable. This is because the 
subtracted backgrounds obtained with the average raw 
signal multiplicity (pi) are too large. In fact, the cumu- 
lant with pure flow correlation and jet-flow cross-terms 
removed is very similar to that shown in Fig. [DJc) for the 
uniform background case. 



IV. SUMMARY 

We have described two analysis methods for three- 
particle azimuthal correlations between a high transverse 
momentum trigger particle and two softer associated 
particles: the jet- correlation method and the cumulant 
method. We point out two ways to analyze the cumu- 
lant, the lab- frame cumulant using azimuthal angle in 
the laboratory frame, and the reaction-plane-frame cu- 
mulant using azimuthal angle relative to reaction plane 
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in nucleus-nucleus collisions. We focused on the com- 
parison between the jet-correlation method and the lab- 
frame cumulant method. We studied their differences an- 
alytically for two cases: one with a uniform background 
in azimuth and the other with a background including 
anisotropic flow. The major conclusions from our study 
are as follows: 

(1) Cumulant is a mathematical construction; different 
mathematical representations of the same physics 
origin yield different cumulants. Extreme care 
should be taken in making physics interpretations 
of cumulant results. 

(2) The reaction-plane-frame cumulant method and 
the jet-correlation method are identical if the back- 
ground in the latter is obtained from mixed-events 
without normalization scaling. Since the azimuthal 
angle relative to reaction plane is a natural vari- 
able to use, we advocate the reaction-plane-framc 
method for future angular cumulant analyses in 
nucleus-nucleus collisions. 

(3) The jet-correlation method has the jet-model built- 
in and is designed to study jet-like correlations; its 
interpretation is straightforward. The lab-frame 
cumulant method does not have a particular corre- 
lation model built-in and studies any kind of corre- 
lations; its interpretation is, however, difficult and 
has to involve a physical correlation model. 

(4) For pure anisotropic flow without jet correlation, 
the jet-correlation method yields zero jet signal. 
The lab-frame cumulant method results in irre- 



ducible flow correlation terms and residual flow cor- 
relation due to non-Poisson statistics. 

(5) With both anisotropic flow and jet correlation 
present, the jet-correlation method can construct 
the jet-correlation signal. However, the signal 
is distorted if incorrect background level is sub- 
tracted, due to any mismatch in anisotropies of jet- 
correlated particles and medium background parti- 
cles. The lab-frame cumulant method results in 
jet-flow cross terms, in addition to the irreducible 
flow-correlation terms in (4) and the jet-correlation 
signal. The resultant lab-frame cumulant is very 
complex. 

(6) For narrow jet peaks, both methods can identify 
Mach-cone structures, although the magnitudes 
and shapes of the correlation functions are differ- 
ent. For broad and more or less flat two-particle 
jet-correlations on the away side, as measured in 
RHIC experiments, the jet-correlation method can 
still identify Mach-cone structures, while the cumu- 
lant method fails. 
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